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If microlensing of stars by dark matter has been detected (Alcock et al. 
1993; Aubourg et al. 1993; Udalski et al. 1993; Alcock et al. 1994; Udalski et al. 
1994a,b), then the way is open for the development of new methods in galactic 
astronomy. This series of papers investigates what microlensing can teach us 
about the structure and shape of the dark halo. In this paper we present formulas 
for the microlensing rate, optical depth and event duration distributions for a 
simple set of axisymmetric disk-halo models. The halos are based on the "power- 
law models" (Evans 1993, 1994) which have simple velocity distributions. 

Using these models, we show that there is a large uncertainty in the predicted 
microlensing rate because of uncertainty in the halo parameters. For example, 
models which reproduce the measured galactic observables to within their errors 
still differ in microlensing rate towards the Magellanic Clouds by more than 
a factor of ten. We find that while the more easily computed optical depth 
correlates well with microlensing rate, the ratio of optical depth to rate can vary 
by a factor of two (or greater if the disk is maximal) . Comparison of microlensing 
rates towards the Large and Small Magellanic Clouds (LMC and SMC) and 
M31 can be used to aid determinations of the halo flattening and rotation curve 



slope. For example, the ratio of microlensing rates towards the LMC and SMC is 
~ 0.7-0.8 for E0 halos and ~ 1.0-1.2 for E7 halos (c.f. Sackett & Gould 1993). 
Once the flattening has been established, the ratio of microlensing rates towards 
M31 and the LMC may help to distinguish between models with rising, flat or 
falling rotation curves. Comparison of rates along LMC and galactic bulge lines- 
of-sight gives useful information on the halo core radius, although this may not be 
so easy to extract in practice. Maximal disk models provide substantially smaller 
halo optical depths, shorter event durations and even larger model uncertainties. 

Subject headings: dark matter - Galaxy: structure - gravitational lensing 

1. Introduction 

The recent detection of possible gravitational microlensing events (Alcock et 
al. 1993; Aubourg et al. 1993; Udalski et al. 1993, Alcock et al. 1994; Udalski et al. 
1994a,b) gives hope that at least part of the dark matter content of our galaxy is 
directly accessible to observation. The dark halo, whose extent has been studied 
gravitationally for many years via velocities of stars, gas, and satellites (e.g. Fich 
& Tremaine 1991), contains at least three times (and perhaps more than ten 
times) the mass of the luminous galaxy. Its identity is one of the major unsolved 
problems in astronomy (e.g. Ashman 1992; Primack, Seckel, & Sadoulet 1988). 
While it is possible that the halo consists mostly of exotic non-baryonic elemen- 
tary particles, the idea of Paczyhski (1986) of searching for Massive Compact 
Halo Objects (Machos) in the range 1O -8 M0 to 1O 3 M0 by monitoring millions 
of stars in the LMC may have borne fruit in the experimental programs. 

If the Milky Way halo contains large numbers of Machos, then the gravi- 
tational microlensing experiments now under way should have the potential to 
determine the number and distribution of Machos in the halo, as well as their 
mass distribution. The next step for the microlensing experiments will be to 
gather more events and then translate the number and duration of those events 
into an estimate of the mass fraction / of the dark halo which consists of Machos 
in the relevant mass range. To accomplish this goal a model of the dark halo 
is necessary. In the past simple spherical models with flat rotation curves have 
been considered (Paczynski 1986; Griest 1991; DeRujula, Jetzer & Masso 1991; 
Nemiroff 1991). These have been valuable in estimating the order-of-magnitude 
effects but suffer from at least three important deficiencies: 

(1) The halo may not be spherical. N-body simulations of gravitational collapse 
of collisionless dark matter generically produce axisymmetric or triaxial halos 
(Quinn, et al. 1992; Dubinski & Carlberg 1991; Katz 1991). The recent papers 
of Sackett & Gould (1993) and Frieman & Scoccimarro (1994) have made an 
important start on the study of microlensing effects in flattened halos (see also 
the early work of Jetzer 1991). 

(2) The effect of the galactic disk is ignored. The disk makes a significant con- 
tribution to the local circular speed. Modeling without proper allowance for this 
effect leads us to over-estimate the gravity field and hence the mass of the dark 
halo. That is, since the amount of material in the galactic halo is set by the local 
circular speed, a larger contribution to this speed by the disk, means a smaller 
halo is needed to explain the total circular velocity. 

(3) It is only a simplified view of the data that permits one to regard the rotation 
curve of the Milky Way as flat. In fact, even the sign of the local gradient of the 
rotation law at the sun is not known - Fich, Blitz & Stark (1989) estimate that 



it may be rising or falling by about 30 km/sec outwards from the solar circle to 
a galacto centric radius of 17 kpc. 

In this paper, we make a start towards quantifying and remedying these defects 
by calculating the microlensing rate and optical depth in a set of simple, flexible 
and realistic halo models. We take the power-law models (Evans 1993, 1994, 
hereafter E93, E94) - for which simple and self-consistent distribution functions 
are known - and provide an approximate method to allow for the influence of the 
galactic disk. This enables calculation not just of the optical depth, but also the 
event rate and the distribution of event durations. In this paper, the emphasis is 
on the uncertainties in microlensing predictions and what can be done to reduce 
them. Since there are many possible models of the dark halo consistent with 
current observations, there is a substantial scatter in the predicted microlensing 
rate - and this in turn contributes to an uncertainty in the measurement of the 
fraction of the halo consisting of Machos. Of course, this is directly relevant to 
the difficult but important question of whether the observed microlensing rates 
can rule out or support the existence of non-baryonic matter in the halo. 

We show that the often calculated optical depth is a useful predictor of the 
microlensing rate to a factor of ~ 2. For more accurate work and to predict event 
durations, a distribution of Macho velocities is needed. By considering a large 
range of model parameters consistent with observations, we find the microlensing 
optical depth can vary by a factor of six or more. If the disk of the Milky 
Way is "maximal" - in the sense that it provides almost all of the Galacto- 
centric acceleration - then a much smaller halo is required. This gives an even 
larger spread in predicted optical depth, rate, and event durations. An attractive 
possibility is to use the measured microlensing rates toward different sources (e.g. 
LMC, SMC, M31 and the galactic bulge) to determine some of the halo and disk 
parameters, thereby providing a new tool for the study of galactic structure and 
at the same time reducing the halo uncertainty in the measurement of /. We 
corroborate the Sackett & Gould (1993) prediction that the ratio of LMC to 
SMC optical depth is a robust indicator of the flattening of the dark halo - and 
extend it by showing that the ratio of microlensing rates distinguishes flatness as 
well. We predict that the ratio of rates toward M31 and the LMC may be enable 
us to discover whether the rotation curve is rising or falling (or equivalently the 
extent of the dark halo). A comparison of the bulge and LMC rates can provide 
information on the halo core radius. 

The plan of the paper is as follows: In §2 we review the axisymmetric models 
and present formulas for the optical depth, microlensing rate, and event dura- 
tions. In §3 we discuss the halo parameters and their allowed range and explain 
how to take into account the effect of the galactic disk. In §4 we compare op- 
tical depth and rate, and discuss the uncertainties in microlensing rates due to 
uncertainties in the parameters. In §5 we discuss reducing those uncertainties by 
comparing results along different lines-of-sight, in §6 we discuss distributions of 
event durations, and in 57 we summarize our conclusions. 
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2. Axisymmetric models 

The primary goal of the galactic gravitational microlensing experiments is 
to determine the mass of the dark halo in Machos. The experiments search for 
Machos by monitoring millions of stars nightly in the Large Magellanic Cloud 
(LMC), the Small Magellanic Cloud (SMC), the galactic bulge, and perhaps in the 
future in the M31 galaxy (Crotts 1992; Baillon et al. 1993). If the dark halo con- 
tains large numbers of Machos, occasionally one passes close to the observer-star 
line-of-sight and acts as a gravitational lens, causing a time-dependent magnifi- 
cation of the stellar image. The resulting lightcurve is determined by only a few 



quantities such as the distance s from us to the Macho, Macho mass m, Macho 
transverse velocity v±, and impact parameter b. The magnification A(t) as a 
function of time t is given by 

A(t) = (u 2 + 2)/[u(u 2 + 4) 1 /% 

u(t) = b/R e = (u 2 min + u 2 (t - to) 2 ) 1 / 2 , (1) 

R e = (2/c)[Gms(l - s/L)] 1 / 2 , 

where the peak magnification A m&x is given by inverting u m [ n = u(A max ), uj = 
v±/R e , and L is the distance to the star. Experimentally, microlensing events 
are characterized by the maximum magnification A max , the time of the peak to, 
and the duration of the event t, where t = 2/oj. Also used in the literature as an 
"event duration" is t e = t(urn — uj^) 1 ' 2 . This is time for which A > Aj>, with 
Ay = A(ut)- Using A^ = 1.34 corresponds to ut — 1 which is the time inside 
the Einstein radius R e . A more thorough discussion can be found in in Paczyhski 
(1986) and Griest (1991). 

An observing team measures the number and duration of microlensing events. 
The number of observed events is proportional to the number of stars monitored, 
the duration of the experiment, the experimental efficiency, and the rate at which 
microlensing occurs. The primary observables are the optical depth r, the rate 
T and the average duration of the events (t). They are related by 



r = r(t e ) = ^r<t). (2) 

If the distribution of Machos masses n(m) were a delta function, then F would 

be T = Ti(77i/Mq)^ 1 ' 2 1 where Ti is the rate with m = M Q and is independent 
of mass. For a general normalized mass distribution F = r] m ri, where the mass 
integral is 

, m =f dm „ im) ( m /M r^ (3) 



This also implies that (t) = (t) 1 rj, 
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The optical depth is the number of Machos inside the microlensing tube with 
radius UTR e (s) and length L. It depends only on the density of Machos p 



t = I P^-TxulRlis^ds. (4) 

m 



Unlike the rate or average duration, it is independent of the Macho mass distri- 
bution. For this reason it also contains no direct information about the mass of 
the lensed objects. Note also that since detection efficiencies depend upon the 
duration of events, it is important to have models which predict durations. Now 
to find the rate at which Machos enter the microlensing tube requires knowledge 



of the distribution of velocities all along the tube. So the rate and the distri- 
bution of event durations are hard to calculate because they require the entire 
phase space distribution function (DF) F(y, x). The differential rate is given by: 



1 o 

dF = — F(v,x) cos 9uTR e v±d vdxda, (5) 

m 



where the angles and notation are defined in Griest (1991). 

What makes the calculation particularly difficult is that the DF cannot be 
prescribed arbitrarily as a Maxwellian, for instance. This is because the Machos 
are collisionless, so the DF is constrained to obey the collisionless Boltzmann 
equation. By Jeans' theorem, this implies that the DF depends only on the iso- 
lating integrals of motion (see Binney & Tremaine 1987, p. 220). Self-consistent 
solutions for distributions of velocities that build flattened halo models are scarce. 
The largest known set of axisymmetric models with simple DFs are the "power- 
law galaxies" (E93, E94). These form the basis for the exploration of microlensing 
in this paper, allowing us to go beyond simple spherical models. Note, however, 
that all these models are axisymmetric and oblate, while N-body simulations 
suggest that halos may well be triaxial. Exploration of triaxial models will be 
done in a future paper. 

The parameters of the power-law models are: 

(1) The core radius R c , which measures the scale at which the density law begins 
to soften. 

(2) The flattening parameter q, which is the axis ratio of the concentric equipo- 
tential spheroids, with q = 1 representing a spherical (E0) halo and q ~ 0.7 
representing an ellipticity of about E6. The "isophotal" ellipticity of the dark 
halo is a function of q, as well as other parameters of the model [see E94, eq. 
(2.9)]. 

(3) The parameter (3, which determines whether the rotation curve asymptotically 
rises, falls or is flat. At large distances R in the equatorial plane, the rotation 

velocity u c i rc ~ R~@. So f3 = corresponds to a flat rotation curve, while f3 < 
is a rising rotation curve and (3 > is falling. 

(4)The solar radius Rq, which is the distance of the Sun from the galactic center. 

(5) Finally, the normalization velocity vq, which determines the overall depth 
of the potential well and hence the typical velocities of Machos in the halo. In 
the limit (3 — 0, q — 1 and large R (spherical halo with a flat rotation curve) 

V0 = fcirc- 

Using z as the height above the equatorial plane, the potential of the power- 
law models is 



v 2 R p c /P 



(R 2 c + R 2 + z 2 q- 2 )P/ r 



*= < 



(6) 
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\og{R z c + R z + z z q- z ), if = 0, 



and the mass density is 

v$Rg R 2 C (1 + 2g 2 ) + R 2 (l - (3q 2 ) + ^ 2 (2 - (1 + /j)q- 2 ) 
P 4nGq 2 (R 2 + R 2 + Z 2 q- 2 )C 3 +^/ 2 • ^ J 

The DF corresponding to this potential-density pair is 

' AL 2 |£| 4 //3-3/2 + S| S |4//J-l/2 + C | E |2//3-l/2 ; if ^ _^ 0? 

F(E,L Z ) = < 

k AL 2 exp(4£/^) + Bexp(4£?/u§) + Cexp(2£/i; 2 ), if /? = 0. 

(8) 
where the constants A, 5, and C are given in E93 and E94. As required by Jeans' 
theorem, the DFs depend only on the isolating integrals of motion, namely the 
relative energy per unit mass E = \1/ — ^i> 2 , and the angular momentum per unit 
mass about the symmetry axis L z . The circular velocity in the equatorial plane 
is 

v 2 . = 4RcR 2 (g) 

CirC ( J R2 + J R 2 )(/3+2)/2- ^ 

Note that the limit q = 1, (3 — and R c = recovers the standard singular 
isothermal sphere used by Paczyhski. Allowing a core radius gives 



vfi R 2 + 3RJ 
4nG(R 2 + R 2 ) 2 ' 



(10) 



This differs from the cored isothermal sphere considered by Griest (1991) in sev- 
eral ways. First, the rotation curve approaches its asymptotic value more quickly. 
Second, the DF given by the q = 1, (3 — limit of equation (8) is self-consistent, 
whereas Griest (1991) assumed an approximate Maxwellian distribution of ve- 
locities. 

So far we have only modeled the dark halo. However, in the standard model, 
a substantial fraction (~ 40%) of the centripetal force at the solar radius derives 
from the disk stars. This is represented by a thin exponential disk with a scale 
length of Rd = 3.5 kpc, normalized to a surface density of Eo = 50M Q pc~ 2 at the 
solar radius (Gilmore et al. 1989; Gould 1990). It is possible that the disk of our 
galaxy is substantially larger than the canonical value. (Oort 1960; Bahcall 1984; 
Kuijken and Gilmore 1989; Gould 1990). Recent microlensing results (Alcock, 
et al. 1994; Udalski, et al. 1994) as well as studies of the optical rotation curves 
of external galaxies (Buchhorn 1992; Kent 1992) may suggest this. We consider 
such a "maximal disk" by taking Eo = lOOM0pc~ 2 . The rotation velocity added 
in quadrature is thus (Freeman 1970; Binney & Tremaine 1987, p. 77) 

v 2 disk = A-KG^hy 2 [I (y)K (y) - h{y)Ki{y)\ , (11) 

where y = R/(2R,i), and the I n and K n are modified Bessel functions. Note 
that in adding a contribution from the disk to the local circular velocity, we have 



sacrificed self-consistency. Really, we should find the DF of the power-law halo 
in the combined potential field of both disk and halo - instead, we use the DF 
(8). As has been argued elsewhere (Evans & Jijina, 1994), this is an reasonable 
approximation for the LMC, SMC and M31, where microlensing typically occurs 
at heights above the equatorial plane of many kpc. 

In this paper, our aim is to estimate the contribution of the Galactic halo to 
microlensing. Of course, this is not the only possible source of deflectors. Towards 
the LMC, there is the possibility of microlensing by the LMC dark halo or disk 
(Gould 1993b, Sahu 1994). The optical depth is ~ 2.5 x 10~ 7 for microlensing by 

LMC halo lenses and ~ 0.09 x 10~ 7 for LMC disk lenses. The Galactic halo makes 
a contribution that is roughly three times greater and so is the dominant source 
of lenses. However, this is not the case for lines of sight towards M31. The optical 
depth is dominated by Machos in the halo and disk of M31 (Crotts 1992, Gould 
1993a). Crotts (1992) estimates that the halo of our own Galaxy contributes just 
20% to the total optical depth. Microlensing towards the Galactic bulge poses 
perhaps the hardest problems of separating the contributions of different deflector 
populations. Bulge stars can undergo microlensing not only by halo Machos, but 
also by other bulge and disk stars (Griest et al. 1991, Paczyhski 1991, Kiraga 
& Paczyhski 1994). At Baade's Window, the optical depth is ~ 6.3 x 10~ 7 for 

microlensing by bulge lenses, ~ 5.0 x 10~ 7 for disk lenses. The dark halo only 
makes an important contribution if the core radius is small. 

We are now in a position to calculate the microlensing observables - the 
optical depth, rate and average duration of events. They can be found using 
equations 2, 3, 4 and 7. The results are single quadratures and readily evaluated 
on the computer. They are displayed in Appendix A. In Appendix B, we give 

the differential microlensing rate dT/dt, where t is defined just after equation (1). 

The probability of obtaining an event of duration t is just (dT/dt)/T. 

3. Range of models 

In order to explore the scatter in microlensing observables, we build a set 
of halo models which span the observationally allowed range. The power-law 
galaxy models allow us to vary the flattening, core radius and rotation law, and 
we consider both canonical and maximal disks. For each parameter in the model, 
we therefore find the range permitted by the observations. Then, several values 
of each parameter are chosen to represent the range. We also ensure that each set 
of parameters gives a model consistent with the measured Milky Way rotation 
curve. So, we study the statistical properties of an ensemble of models, each one 
of which is a plausible representation of the dark halo of the Milky Way. 

For the dark halo flattening, little is known. So the entire range of flatten- 
ing allowed by the power-law models is examined. This varies between E0 or 
spherical (q = 1) and roughly E6 or E7 (depending on j3). The core radius of the 
dark halo is also uncertain - Bahcall, Schmidt & Soneria (1983) estimate R c as 
2 kpc from star count data, while Caldwell & Ostriker (1981) suggest 10 kpc. If 
the disk is maximal, values as large as 20 kpc are possible. We consider values 
of 2 kpc, 5 kpc, 10 kpc, and 20 kpc. The parameter (3 determines the slope of 
asymptotic circular velocity. Between Rq and 2Rq, the circular velocity is prob- 
ably within 10 — 15% of the I.A.U value of 220 km/s, but whether the measured 
HI rotational velocities rise or fall with R depends upon estimates of the solar 
position Rq and the local circular speed v c irc(Ro) (see Fich et al. 1989; Jones et 
al. 1993). Beyond 20 kpc, little is known directly, though arguments based on 
the kinematics of distant satellite galaxies support the idea of a relatively flat 



rotation curve out to an unknown cut-off (Fich & Tremaine 1991). However, 
current theories of galaxy formation tend to favor the alternative view that dark 
halos extend indefinitely, fading into structure on larger scales. So, we do not 
consider a halo cut-off in this paper - it would add yet another poorly known 
parameter to our model. We investigate power-law halos with j3 = -0.2, 0, and 
0.2. These correspond to rotation curves which rise by ~15%, are flat, or fall 
by ~15% between the solar radius and twice the solar radius, depending a little 
upon R c . 

The value of the solar radius Rq has been reviewed by Reid (1989). He 
shows that most recent determinations lie between 7 kpc and 9 kpc, with 7.7 
kpc being his preferred value. This differs considerably from the IAU value of 
8.5 kpc (Kerr & Lynden-Bell 1986). We examine the values _Rn = 7, 8, and 9 
kpc. Finally, perhaps the single most important parameter is the normalization 
velocity vq. Given our fixed disk contribution to the total rotation law, the 
parameter vq is now specified once we settle upon a choice for w c i rc (i?o)- Merrifield 
(1992) estimates v circ (R ) = 200 ± 10 km/s, Fich, Blitz, & Stark (1989) give 
v c i IC (Ro) = 220 ± 30 km/s, while Rohlfs et al. (1986) give values between 170 
km/s and 200 km/s between Rq = 6 kpc and Ro = 16 kpc. For our ensemble of 
models, we impose the constraint that the total circular velocity lies between 180 
km/s and 250 km/s at Ro and 2Rq. Note that the IAU value is 220 km/s (Kerr 
& Lynden-Bell 1986). We also investigated a more restricted ensemble of models 
with 190 < fdrc(-Ro) < 230 km/s. We find all our results also hold for this more 
restricted ensemble. 

4. Uncertainties in the Rates 

First, let us consider the difference caused by using the optical depth instead 
of the microlensing event rate. The optical depth to microlensing is the mean 
number of Machos in the microlensing tube; that is the number of microlensing 
events taking place at a given moment. It is easy to calculate since it is inde- 
pendent of lens mass and velocity, and only requires knowledge of the density 
distribution p(x). For this reason, it is the most widely estimated quantity. But 
how well does it trace the microlensing rate? 

We are able to answer this question since both the rate V (equation Al) and 
the optical depth r (equation A6) are known for the power-law models. One way 
to test this is to plot (t), which is the ratio of optical depth r and T, (t) = ^t/T, 
for many different models. The average duration (t) is a constant if r and T are 

well-correlated. In Fig. 1, we show histograms of (t) for microlensing towards 
the LMC, SMC and M31 for our ensemble of models. Figs, la-c demonstrate 
that (t) tends to vary by more than a factor of two between models. Figs, ld-f 
show an even larger spread for maximal disk models. Figs. 2a-f show this another 
way by plotting the rate vs the optical depth for the set of models. These plots 
show that (t) is indeed much less model dependent than either r or T. While the 
rate and r vary by more than a factor of ten in these plots, their ratio varies only 
~ 2 for a canonical disk. In fact, we note that the line T oc r 3 ' 2 is a fairly good 
fit to all the models we have considered. Thus the large scatter in (T) seen in 
the maximal disk histograms is mostly just due to the large scatter in rate. (The 



* To the extent that the relation T = arj m f 1 / 2 r 3 / 2 holds, where a is a constant from theory, 

we have that a = / _1/ ' 2 r7~ 1 (4/7r) 3 / 2 (t) r~ 1//2 is independent of the model parameters. 
Thus, if the macho fraction / were known, one could extract the mass integral r\ m from 



rate varies more than the optical depth.) In all the plots we use m = M Q , but for 
an arbitrary mass distribution just scale V by eta m and (t) by r^ 1 . Keep in mind 
that a given experiment can produce only one point in the T, r plane and that 
the primary use of a measurement will be to find / , the Macho fraction. We see 
that for approximate work, the optical depth does a reasonable job of predicting 
the rate. But for more detailed work, especially when efficiencies are involved, 
the difference between rate and optical depth should be kept in mind. We also 
note that the predicted distribution of event durations is found as a differential 
rate (Appendix B). 

Next let us turn to scatter in the predicted microlensing rate caused by un- 
certainties in the halo parameters. Fig. 2 shows that for all lines of sights, there 
is a scatter in the rate of more than a factor of ten for a canonical disk. For 
the LMC, the models with the smallest rate have spherical halos with small core 
radii, falling rotation curves, and small values of i>o, while the models with the 
largest rates have either spherical or flattened halos, but large core radii, rising 
rotation curves, and large values of vq. This is as expected, since any model 
which puts more mass at a large distance in the direction of the LMC will have 
a larger microlensing rate, and a larger optical depth. This is shown in Fig. 3 
in which we plot the optical depth against the rotation velocity at r = 50 kpc. 
The correlation between r and v\{r = 50 kpc), while not perfect, is quite good. 
Note that the mean value of the rotation velocity at Ro is nearly independent of 
the microlensing rate. Figs. 2d-f and Fig. 3c-d show the case of a maximal disk. 
Here we see that the rate and optical depth can be considerably smaller than 
for a canonical disk. Also there is a variation between models of several orders 
of magnitude. This is as expected since in these models the disk is the main 
contributor to the rotation curve at the solar distance. Thus a smaller enclosed 
halo mass is required to match observations, and the halo parameters are poorly 
constrained. 

The halo may only consist of a fraction / of baryonic matter in the form 
of Machos. Thus, a factor of more than ten uncertainty in the predicted rate 
caused by the poorly determined halo parameters makes it difficult to determine 
the allowed amount of non-baryonic dark matter. It is clearly essential to reduce 
the uncertainty. 

5. Reducing Model Uncertainties 

The primary way of reducing the model uncertainties in the microlensing 
observables is to determine the halo parameters. Even within the restricted 
framework of the power-law galaxy models, if f3, vq, q, Rq and R c are known, 
there is still uncertainty in the rate. This is because the DFs equation (8) are 
the simplest consistent with the potential and the density, but are certainly not 
unique. There are still further DFs that depend on non-classical third integrals 
of motion and generate anisotropic velocity distributions. Note, too, that even 
though our models give a plausible representation of the Milky Way, there cer- 
tainly exist other alternatives (see e.g., Frieman & Scoccimarro 1994; Gates & 
Turner 1993; Giudice, Mollerach & Roulet 1994) with different lensing properties. 
And of course, the size of the disk plays a crucial role. 



observables rj m w f 1 / 2 a~ 1 ('i/n) 3 / 2 N e ffE 1 / 2 ('%2 ti/ti)~ 3 ^ 2 , where E is the total exposure, 
N e ff = X) e r ' e i * s ^ ne efficiency at which events of duration i, are recovered, and the sums 
go from 1 to the number of observed microlensing events. For LMC microlensing in our set 
of models we find a s=a 3850 ± 260yr~ 1 . The physical basis for this relationship may simply 
be that the optical depth is proportional to the mass along the line of sight ex v 2 , and the 
rate is proportional to the optical depth times v c . 



One obvious way to determine halo parameters is to use conventional astro- 
nomical techniques - observations of stars, gas and satellites - to fix the solar 
radius and circular speed more accurately. For example, fixing the solar radius 
at 8 kpc, and demanding v c [ rc = 220 km/s ±5% between 8 and 16 kpc reduces 
the spread in microlensing rates toward the LMC from more than a factor of ten 
to a little more than a factor of two (for the canonical disk). Uncertainties in r 

and (t) are reduced similarly. A better determination of the halo core radius by 
stellar observations would also be important. 

However, it is also possible to use the microlensing experiments themselves to 
determine the halo parameters and reduce the model uncertainty. The basic idea 
is to exploit the fact that there are at least four viable lines-of-sight out of the 
Milky Way in which to measure the microlensing rate and average event duration. 
Each line-of-sight (LMC, SMC, M31 and the bulge) offers a different "pencil 
beam" through the dark halo, and so by comparing the rates, optical depths, 
and average durations among the different lines-of-sight information concerning 
the halo shape can be gained. Several of the parameters, such as flattening q 
and asymptotic slope of the rotation law /?, may best be determined this way. 
So microlensing gives us a new probe of the density and velocity structure of the 
dark halo. This is in addition to information on the size of the disk gained via 
microlensing. 

For instance, a scatter plot of the ratio of LMC and SMC rates vs the LMC 
and SMC average durations is shown in Fig. 4a. The models clearly fall into 
two distinct groups. Those models marked with a circle all have round halos 
(E0), while those with a square are flattened to roughly E6. Thus the ratio of 
LMC rate to SMC rate is a excellent indicator of halo flattening. This effect was 
first discovered - using optical depth rather than microlensing rate - by Sackett 
& Gould (1993). Frieman & Scoccimarro (1994) have recently cautioned that 
the robustness of this diagnostic may be lost if the halo is tilted with respect to 
the disc - although such a configuration cannot be a long-lasting equilibrium. 
So, the halo flattening can probably be determined if enough events are found 
to allow accurate measurement of the SMC microlensing rate. Figs. 4b and 4c 
show the rate ratio for M31/LMC, and M31/SMC. While separation of flattened 
models is still evident, one sees from the figures that it is the LMC and SMC 
position relative to the halo axis of symmetry that make the measurement of 
the flattening so easy. Note again, that in an experiment one measures only one 
LMC rate (and optical depth) and one SMC rate, and so gets only one point in 
any of these scatter plots. It is also interesting to observe from Fig. 4a that the 
model uncertainties in the LMC/SMC rate ratio are much greater for flattened 
halos than for spherical halos. The case of a maximal disk is not shown, since it 
looks almost identical to Fig. 4. 

Can we use microlensing to determine whether the halo has a rising or falling 
rotation curve? The LMC and SMC are at nearly the same distances (50 and 60 
kpc), so it is natural to expect the ratio of M31 to LMC microlensing to be the 
most useful discriminant. Note that rate ratios are convenient to use, because the 
magnitude of any rate always contains the unknown parameter /. In Fig. 5 we 
plot the M31/LMC rate ratio vs the LMC rate for the set of models above, with 
triangles for (3 = 0.2 (falling rotation curve), circles for (3 — (asymptotically 
flat rotation curves), and stars for (3 = —0.2 (rising rotation curve). In Fig. 5a, 
all models are plotted, while in Fig. 5b and Fig. 5c only models with spherical 
(q = 1) and flattened (q = 0.71 or q = 0.78) halos respectively are shown. 
In Fig. 5a some separation of models with different values of (3 is evident but 
there is substantial ambiguity, which would make a direct estimate of (3 using 
this method difficult. However, suppose that we have already determined the 
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halo flattening by use of the ratio Tlmc/^smc- Then, as shown in Figs. 5b 
and 5c for a canonical disk, a fairly clear separation of rising, falling, and flat 
rotation curve parameter can be accomplished. Thus, the ambiguity seen in 
Fig. 5a is largely removed when models with different flattenings are plotted 
separately. The exception is some overlap between models with R c = 2 kpc and 
R c = 20 kpc and different values of (3. This ambiguity is probably removable 
as discussed below. The case of a maximal disk is not displayed, as it is very 
similar. So, the asymptotic form of the rotation law, or equivalently (3, can 
probably be determined from the M31/LMC rate ratio once q is known. Keep in 
mind, however, the caveats mentioned in §3 concerning our M31 rate calculation, 
which may result in corrections which modify this effect. If halo microlensing 
can be distinguished from M31 microlensing, a measurement of (3 should then be 
possible. 

Next, can we determine the halo core radius _R C ? The parameter R c affects 
mainly the inner portion of the halo and overall normalization of the halo mass. 
This overall normalization is mixed in with vo and /, and so the best hope 
in determining R c is probably a comparison of the bulge with a more distance 
source such as the LMC. Here we have the problems mentioned in §3 concerning 
bulge microlensing; our modelling of the distribution of velocities is not adequate 
along the disk. But, the optical depth is independent of the velocities and will 
give some indication of the rate. Even so, our calculations do not give the total 
optical depth towards the bulge, merely the contribution of the optical depth 
from the halo. 

In Fig. 6, we plot the LMC/bulge optical depth ratio vs the bulge optical 
depth, where triangles indicate R c = 2 kpc, boxes indicate R c = 5 kpc, and 
stars indicate R c = 10 kpc. A reasonably clean separation is obtained when 
this ratio is plotted for all the models (not shown). In Fig. 6, this separation 
is made clear-cut, if one supposes (3 and q have already been measured by the 
methods above. The R c = 20 kpc models have an LMC/bulge ratio of greater 
than 10, and are very easily distinguished even with no prior knowledge of (3 
and q. (They fall off the top of the plots in Fig. 6). Even if (3 and q are not 
known, the separation is quite good if the value of the solar radius R$ is held 
fixed. So, a better determination of Rq by non-microlensing means can allow a 
clearer separation of the effect of the halo core radius. The case of a maximal 
disk is not shown since it gives very similar results. 

6. Distribution of Event Durations 

Since the duration of a microlensing event is proportional to the Einstein 
radius (oc to 1 ' 2 ), the duration of an event gives information about the mass of lens 
which caused it. In trying to understand the nature of the objects responsible for 
the observed microlensing, this is important information. But the duration also 
depends upon the unknown lens velocity and distance. Thus, a given mass Macho 
can cause a wide distribution of event durations. This distribution must be used 
statistically to infer probable masses from observed durations. Using the DF's 
(equation (8)), the distribution of event durations can be found. The formula and 

definitions are given in Appendix B. In Fig. 7, we show several t distributions. 
One sees that different halo parameters give quite different distributions. It is the 
average of these distributions (?) that is shown in the histograms in Fig. 1. Fig. 7 
shows that, as expected, uncertainty in the halo model will lead to additional 
uncertainty in determining the masses of the lensing objects. The curves labeled 
(a), (b), and (c) are canonical disk cases with various choices of halo parameters, 
while curve (d) shows a maximal disk example. We also note that the scaling 
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introduced in Griest (1991) works fairly well for models we considered. That 

is, by scaling the t axis by (t) , and the dT/dt axis by (t), all the curves are 
found to lie roughly on top of each other. This means that for power law galaxy 
models along a given line-of-sight, the shape of the distribution is much more 
model independent than peak value. 

In a future paper we plan to explore further the information that can be 
extracted from event duration distributions, and include other possibilities such 
triaxiality, streaming motion, etc. 

7. Conclusions 

This paper has shown how to exploit the power-law galaxy models (E93, E94) 
as simple, flexible and realistic representations of the dark halo. These models 
have the advantage of simple and analytic phase space distribution functions and 
therefore permit accurate calculation of the optical depth, microlensing rate and 
average event duration. We provide formulae for these quantities as a function of 
the halo parameters and source distance and direction (Appendix A). The distri- 
bution of event timescales is presented in Appendix B. We apply our formulae to 
study microlensing towards the Large and Small Magellanic Clouds (LMC and 
SMC), the galactic bulge, and the M31 disk galaxy. We find that: 

(1) For a canonical disk, the optical depth is a reasonable indicator of the mi- 
crolensing rate to within a factor of two. This is important, because the optical 
depth is much easier to calculate than the rate and probably will continue to be 
widely used by investigators. For more accurate work, as well as for derivations 
of the distribution of durations, galaxy modeling with distribution functions is 
crucial. For a maximal disk the agreement between optical depth and rate is less 
robust, though the relation T oc m ' r ' seems to hold. 

(2) The evaluation of the fraction / of the halo consisting of Machos is hampered 
by the uncertainties in the galactic constants, such as the shape of the rotation 
law and the flattening of the dark halo. For a realistic set of halo models, we 
found rates toward the LMC and SMC can vary by more than a factor of ten 
from model to model for a canonical disk, and by several orders of magnitude for 
a maximal disk. Left unaddressed, this model uncertainty will thwart accurate 
determination of /. 

(3) An attractive way of reducing the uncertainty - which simultaneously opens 
up a new method in galactic astronomy - is to use microlensing to explore the 
shape and structure of the dark halo. This has also been realised by Sackett & 
Gould (1993), who showed that the ratios of optical depth towards the LMC and 
SMC is a robust indicator of the flattening of the dark halo. We confirm this 
result by showing that the ratios of the event rates also distinguish flatness. In 
particular, the ratio of microlensing rates towards the LMC and SMC is ~ 0.7—0.8 
for E0 halos and ~ 1.0 — 1.2 for E7 halos. This is true for both canonical and 
maximal disk models. Once the flattening has been established, the asymptotic 
slope of the rotation curve (3 might be determined using the M31/LMC rate ratio. 
The LMC/bulge ratio contains important information on the halo core radius. 
We caution that this may not be easy to extract, as the dark halo is probably 
not the dominant source of lenses towards the bulge. 

In summary, the discovery of a dark halo consisting of a significant fraction 
of Machos is only the starting point for an exploration of the halo characteristics 
which microlensing can help determine. 
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Appendix A 

In this appendix, we give the formulae for the microlensing rate and optical 
depth for the general flattened halo model described in the text (equations 4-7) . 
The total rate F of microlensing in a power-law halo with model parameters 
(3, vq, R c , Rq, and q is 

_ C u T v$Rf /2 (P + 2)(1 - q 2 ) F{np) 
y/2irM/M Q 2 C ( ? 2 v /^L 1 /2+3/3/2 T(^) 1 



Cqu t vlR 2 c +3f3/2 (f3 + 2)F(n p ) 
2nM/M Q cq 2 V=j3L 5 / 2 +W 2 F(dp) 



h (Al) 

Cqut vlR?'\2 - ffi F{n p - 2/\(3\) 
'2nM/M Q Cy f=PlM*+W* T(d p -2/\p\) 3 ' 

Here, Cq = 1/^/GMq, F(x) is the gamma function, and the integrals I{ are 

f ds y/s(l-s)(A's 2 + B's + C) 

I1= J 



1 



(D's 2 + Es' + F') 2 +W/2 


d,Sy/s(l — s) 


(D's 2 + Es' + F') 2 + 3 ^/ 4 


dSy/s(l — S) 



^2 — / JT^To , ~, , 7-.m94_3«/zl ( A2 ) 



1 

h = 





(D's 2 + Es' + F') 1 +W/*' 



with 



A' =3 cos 6, B' = — 6Ro cosbcosi/L, 
C =2R\jL 2 + Rl cos 2 £/L 2 + R 2 Q sin 2 £ sin 2 b/L 2 , 
D' = cos 2 b + q~ 2 sin 2 6, E' = -2R cos b cos £/L 
F'=(R 2 C + R 2 )/L 2 . 



(A3) 



The quantities 6, / are the galactic coordinates of the source star, L is the source 
distance, G is Newton's constant, and c is the speed of light. The constants np 
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and dp have a different form according to whether (5 is positive or negative 



( -4 3 

T"2- if/3<0 - 



n /3 = \ 



(A4) 



L ^ + 2, if /? > 0, 



da = < 



r -4 

7 



-1, if /? < 0, 



4 5 
U + 2< if " >0 ' 



(A5) 



In the limit /3 — ► (the case of an asymptotically flat rotation curve), the ex- 
pression for the rate follows from the above by systematic use of the formula 
F(x + l/2)/r(x) — > y/x as x — > oo. The optical depth r is 



,2p/3„,2 



r 



itfi#i4 / , s(l-s)(AV + B"s + C")ds 



cV^ 7 (D's 2 + E's + F')(P+*)/ 2 ' 
o 



(A6) 



where 



A" = (1 - /3q 2 ) cos 2 6 + (2 - (1 + /3)g- 2 ) sin 2 6, 
B" = -2(l-(3q 2 )R cosbcos£/L, 
C" = (R 2 c (l + 2q 2 )+R 2 (l-(3q 2 ))/L 2 . 



(A7) 



The quadratures are straightforward to evaluate on the computer. 



Appendix B 



The distribution of event durations is important for finding the mass of 
the lensing objects. It is given by the normalized differential microlensing rate 
(dT/dt)/T, where t = 2R e /v±, and v± is the speed of the Macho perpendicu- 
lar to the line-of-sight. The time the Macho spends inside the Einstein radius, 
t e = (wj, — w 2 nin ) 1 ' 2 t, where u m i n is defined in equation (1), and ut = 1. The 
average duration is related to the average t by (t e ) = j (t). In many cases 

it is advantageous to use distributions in t, since they are independent of the 
amplifications. 
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For the model described in the text, we find: 



dF 

di 



ut ( Is- 



nc 2 \ R^t 4 



((3 + 2) |/3| 1+4 /' 3 (q- 2 - 1) [aiG' Ji + a 2 H' J 2 ] 



+ 



u T ( L 4 \ \f3\ 1+4 /P(f3 + 2) 



ttc 2 \R 2 i 4 



a\J\ 



+ 8 5(sj) l,3|I+2/ " (2 - ? " (1 + Wa3j3 



where, 



(Bl) 



a i 



0-2 



a 3 



-1-4 

1 + 1 
4(1 + 4) 

_1_ 2 



1 + 



/3<0 
/3>0, 



/3<0 
/3>0, 



(B2) 



and the integrals Jj are 



J 2 



dss 2 (l-s) 2 
dss 3 (l-sf 



K' 



H 1 



g 1 (D's 2 + E's + F')f 3 / 2 0i 



f 2 + |W'-(§ 



s(l-a) 



4/0 



*' * 8(1-8) 

g 2 (D's 2 + E's + F')^ 2 g 2 



4//3-1 



J 3 = / dss 2 (l-s)' 



K' 



gs(D's 2 + E's + F')^ 2 03 



— «1 -a 



2/0 



(B3) 



If /3 < 0, the integrals are evaluated over the interval [0, 1]. If (3 > 0, then we 
must restrict the domain of integration by 



p 8/3Ls(l- s) (m/M Q ) (D's 2 + E's + F'fl 2 



{v q cC q ) 2 {R c /LY 



(B4) 



The constants A', B', C, D' and E' are given in Appendix A. The additional 
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constants are 



Rq . ,\ TT i 8 m 



G' = ( —— cos b sin / ) , H' 



L(cC o ) 2 M ' 
2 i /d„\/3 /.. f\ 2 






(B5) 
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Figure Captions 

Figure 1: Histograms of the average duration (tj = —r/Y for the ensemble of 
halo models discussed in the text. Part (a) is for the LMC, (b) is for the SMC, 
and (c) is for M31. If optical depth tracked microlensing rate perfectly each 
histogram would be a delta function. Parts (d)-(f) are the same for a maximal 
disk model. Note all plots are for m = 1M©; scale by rf^ for other masses 
(Equation (3)). 

Figure 2: Scatter plots of microlensing rate vs optical depth for the ensemble 
of models discussed in the text. Part (a) is for the LMC, (b) is for the SMC, and 
(c) is for M31. Each point represent a consistent model of the dark halo. Parts 
(d)-(f) are the same for a maximal disk model. All event rates scale Y oc rj m . 

Figure 3: Scatter plots of optical depth vs i^(50 kpc), the square of the total 
rotation velocity at 50 kpc in the galactic plane. The mass of the Galaxy interior 
to this distance is proportional to this squared velocity. Parts (a) and (b) are for 
a canonical disk, while parts (c) and (d) are for a maximal disk. 

Figure 4: Finding the flattening parameter q. Scatter plots of the ratio of 
rates vs. the ratio of event durations. The circles represent halo models which 
are spherical (q = 1), while the squares represent flattened halos (q = .71 for 
(3 = 0, -0.1; q = 0.78 for (3 = 0.1). part (a) is for LMC/SMC and shows clear 
separation of spherical and flattened halos. Part (b) is M31/LMC, and part (c) 
isM31/SMC. 

Figure 5: Finding the asymptotic slope (3. Scatter plots of the M31 rate divided 
by the LMC rate vs the LMC rate. The stars represent halo models with (3 = —0.2 
(rising rotation curve), the circle models with (3 — (flat), and the triangles 
models with (3 = 0.2 (falling). Part (a) shows all models, while part (b) shows 
only spherical models and part (c) only the flattened models. Separation of the 
models becomes easier if the flattening is known. The line of ambiguity in some 
panels is due to R c = 20 kpc models, which can be distinguished as shown in 
Figure 6. All event rates scale Y oc rj m . 

Figure 6: Finding the core radius R c . Each panel shows models of definite values 
of j3 (rotation curve slope) and q (flattening) . The separation between models is 
quite good. The triangles represent R c = 2 kpc, the squares represent R c = 5 
kpc, and the stars represent R c = 10 kpc. Panels marked E0 are for spherical 
halos, while those marked E6 are for flattened halos. Models with R c = 20 kpc 
were also considered but they are easily distinguished since they typically have 
r(LMC)/r(bul) > 10 and therefore fall off the top of the figures. 

Figure 7: Examples of LMC t distributions for various model parameters. The 
integral under each distribution is unity. The curve marked (a) is for a "standard" 
spherical halo ((3 = 0, q — 1, R c = 5 kpc, vq = 200 km/sec, Rq = 8.5 kpc, and 
r = 1.64 x 10 -6 events/yr). Curve (b) has a shorter average duration ((3 = —0.2, 
q = 1, R c = 5 kpc, v = 200 km/sec, #0 = 8.5 kpc, and T = 3.9 x 10~ 6 
events/yr). Curve (c) has a longer average duration ((3 = 0.2, q = 0.78, R c = 10 
kpc, vq = 210 km/sec, R$ = 8.5 kpc, and Y = 1.24 x 10~ 6 events/yr). Finally 
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curve (d) has a maximal disk, which greatly reduces the amount of halo material 
(ft = 0, q = 1, R c = 20 kpc, v = 90 km/sec, #0 = 7 kpc, and T = 9.37 x 10" 8 
events/yr). The average of each distribution is (t). All event rates scale V oc r\ m . 
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Abstract 



If microlensing of stars by dark matter has been detected (Alcock et al. 
1993; Aubourg et al. 1993; Udalski et al. 1993; Alcock et al. 1994; Udalski et al. 
1994a,b), then the way is open for the development of new methods in galactic 
astronomy. This series of papers investigates what microlensing can teach us 
about the structure and shape of the dark halo. In this paper we present formulas 
for the microlensing rate, optical depth and event duration distributions for a 
simple set of axisymmetric disk-halo models. The halos are based on the "power- 
law models" (Evans 1993, 1994) which have simple velocity distributions. 

Using these models, we show that there is a large uncertainty in the predicted 
microlensing rate because of uncertainty in the halo parameters. For example, 
models which reproduce the measured galactic observables to within their errors 
still differ in microlensing rate towards the Magellanic Clouds by more than 
a factor of ten. We find that while the more easily computed optical depth 
correlates well with microlensing rate, the ratio of optical depth to rate can vary 
by a factor of two (or greater if the disk is maximal). Comparison of microlensing 
rates towards the Large and Small Magellanic Clouds (LMC and SMC) and 
M31 can be used to aid determinations of the halo flattening and rotation curve 



slope. For example, the ratio of microlensing rates towards the LMC and SMC is 
~ 0.7-0.8 for E0 halos and ~ 1.0-1.2 for E7 halos (c.f. Sackett & Gould 1993). 
Once the flattening has been established, the ratio of microlensing rates towards 
M31 and the LMC may help to distinguish between models with rising, flat or 
falling rotation curves. Comparison of rates along LMC and galactic bulge lines- 
of-sight gives useful information on the halo core radius, although this may not be 
so easy to extract in practice. Maximal disk models provide substantially smaller 
halo optical depths, shorter event durations and even larger model uncertainties. 

Subject headings: dark matter - Galaxy: structure - gravitational lensing 

1. Introduction 

The recent detection of possible gravitational microlensing events (Alcock et 
al. 1993; Aubourg et al. 1993; Udalski et al. 1993, Alcock et al. 1994; Udalski et al. 
1994a,b) gives hope that at least part of the dark matter content of our galaxy is 
directly accessible to observation. The dark halo, whose extent has been studied 
gravitationally for many years via velocities of stars, gas, and satellites (e.g. Fich 
& Tremaine 1991), contains at least three times (and perhaps more than ten 
times) the mass of the luminous galaxy. Its identity is one of the major unsolved 
problems in astronomy (e.g. Ashman 1992; Primack, Seckel, & Sadoulet 1988). 
While it is possible that the halo consists mostly of exotic non-baryonic elemen- 
tary particles, the idea of Paczyhski (1986) of searching for Massive Compact 
Halo Objects (Machos) in the range 10~ Mq to 10 Mq by monitoring millions 
of stars in the LMC may have borne fruit in the experimental programs. 

If the Milky Way halo contains large numbers of Machos, then the gravi- 
tational microlensing experiments now under way should have the potential to 
determine the number and distribution of Machos in the halo, as well as their 
mass distribution. The next step for the microlensing experiments will be to 
gather more events and then translate the number and duration of those events 
into an estimate of the mass fraction / of the dark halo which consists of Machos 
in the relevant mass range. To accomplish this goal a model of the dark halo 
is necessary. In the past simple spherical models with flat rotation curves have 
been considered (Paczyhski 1986; Griest 1991; DeRujula, Jetzer & Masso 1991; 
Nemiroff 1991). These have been valuable in estimating the order-of-magnitude 
effects but suffer from at least three important deficiencies: 

(1) The halo may not be spherical. N-body simulations of gravitational collapse 
of collisionless dark matter generically produce axisymmetric or triaxial halos 
(Quinn, et al. 1992; Dubinski & Carlberg 1991; Katz 1991). The recent papers 
of Sackett & Gould (1993) and Frieman & Scoccimarro (1994) have made an 
important start on the study of microlensing effects in flattened halos (see also 
the early work of Jetzer 1991). 

(2) The effect of the galactic disk is ignored. The disk makes a significant con- 
tribution to the local circular speed. Modeling without proper allowance for this 
effect leads us to over-estimate the gravity field and hence the mass of the dark 
halo. That is, since the amount of material in the galactic halo is set by the local 
circular speed, a larger contribution to this speed by the disk, means a smaller 
halo is needed to explain the total circular velocity. 

(3) It is only a simplified view of the data that permits one to regard the rotation 
curve of the Milky Way as flat. In fact, even the sign of the local gradient of the 
rotation law at the sun is not known - Fich, Blitz & Stark (1989) estimate that 



it may be rising or falling by about 30 km/sec outwards from the solar circle to 
a galactocentric radius of 17 kpc. 

In this paper, we make a start towards quantifying and remedying these defects 
by calculating the microlensing rate and optical depth in a set of simple, flexible 
and realistic halo models. We take the power-law models (Evans 1993, 1994, 
hereafter E93, E94) - for which simple and self-consistent distribution functions 
are known - and provide an approximate method to allow for the influence of the 
galactic disk. This enables calculation not just of the optical depth, but also the 
event rate and the distribution of event durations. In this paper, the emphasis is 
on the uncertainties in microlensing predictions and what can be done to reduce 
them. Since there are many possible models of the dark halo consistent with 
current observations, there is a substantial scatter in the predicted microlensing 
rate - and this in turn contributes to an uncertainty in the measurement of the 
fraction of the halo consisting of Machos. Of course, this is directly relevant to 
the difficult but important question of whether the observed microlensing rates 
can rule out or support the existence of non-baryonic matter in the halo. 

We show that the often calculated optical depth is a useful predictor of the 
microlensing rate to a factor of ~ 2. For more accurate work and to predict event 
durations, a distribution of Macho velocities is needed. By considering a large 
range of model parameters consistent with observations, we find the microlensing 
optical depth can vary by a factor of six or more. If the disk of the Milky 
Way is "maximal" - in the sense that it provides almost all of the Galacto- 
centric acceleration - then a much smaller halo is required. This gives an even 
larger spread in predicted optical depth, rate, and event durations. An attractive 
possibility is to use the measured microlensing rates toward different sources (e.g. 
LMC, SMC, M31 and the galactic bulge) to determine some of the halo and disk 
parameters, thereby providing a new tool for the study of galactic structure and 
at the same time reducing the halo uncertainty in the measurement of /. We 
corroborate the Sackett & Gould (1993) prediction that the ratio of LMC to 
SMC optical depth is a robust indicator of the flattening of the dark halo - and 
extend it by showing that the ratio of microlensing rates distinguishes flatness as 
well. We predict that the ratio of rates toward M31 and the LMC may be enable 
us to discover whether the rotation curve is rising or falling (or equivalently the 
extent of the dark halo). A comparison of the bulge and LMC rates can provide 
information on the halo core radius. 

The plan of the paper is as follows: In §2 we review the axisymmetric models 
and present formulas for the optical depth, microlensing rate, and event dura- 
tions. In §3 we discuss the halo parameters and their allowed range and explain 
how to take into account the effect of the galactic disk. In §4 we compare op- 
tical depth and rate, and discuss the uncertainties in microlensing rates due to 
uncertainties in the parameters. In §5 we discuss reducing those uncertainties by 
comparing results along different lines-of-sight, in §6 we discuss distributions of 
event durations, and in §7 we summarize our conclusions. 

2. Axisymmetric models 

The primary goal of the galactic gravitational microlensing experiments is 
to determine the mass of the dark halo in Machos. The experiments search for 
Machos by monitoring millions of stars nightly in the Large Magellanic Cloud 
(LMC), the Small Magellanic Cloud (SMC), the galactic bulge, and perhaps in the 
future in the M31 galaxy (Crotts 1992; Baillon et al. 1993). If the dark halo con- 
tains large numbers of Machos, occasionally one passes close to the observer-star 
line-of-sight and acts as a gravitational lens, causing a time-dependent magnifi- 
cation of the stellar image. The resulting lightcurve is determined by only a few 



quantities such as the distance s from us to the Macho, Macho mass m, Macho 
transverse velocity vj_, and impact parameter b. The magnification A{t) as a 
function of time t is given by 



A(t) = (u 2 + 2)/[u(u 2 + 4) 1 / 2 ] 
i(t) = b/R e = (u 2 min + u 2 (l 
R e = {2/c)[Gms(l-s/L)]^ 2 , 



U(t) = b/Re = « in + U?{t - t ) 2 ) 1/2 , (1) 



where the peak magnification A max is given by inverting u m \ n = u(A max ), u = 
vj_/R e , and L is the distance to the star. Experimentally, microlensing events 
are characterized by the maximum magnification A max , the time of the peak to, 
and the duration of the event t, where t = 2/u. Also used in the literature as an 
"event duration" is t e = t(u T — u^m) • This is time for which A > At, with 
At = A{ut)- Using At = 1.34 corresponds to ut = 1 which is the time inside 
the Einstein radius R e . A more thorough discussion can be found in in Paczyhski 
(1986) and Griest (1991). 

An observing team measures the number and duration of microlensing events. 
The number of observed events is proportional to the number of stars monitored, 
the duration of the experiment, the experimental efficiency, and the rate at which 
microlensing occurs. The primary observables are the optical depth r, the rate 

r and the average duration of the events (tj. They are related by 



r = T(t e ) = -T(t). (2) 

If the distribution of Machos masses n(m) were a delta function, then T would 

be r = T\(m I Mq)~ I , where Ti is the rate with m = Mq and is independent 
of mass. For a general normalized mass distribution Y = rj m Y\, where the mass 
integral is 



r/r, 



/ dmn(m)(m/M Q )~ 1/2 . (3) 



This also implies that ft) = (t\ r] m . 

The optical depth is the number of Machos inside the microlensing tube with 
radius UTR e {s) and length L. It depends only on the density of Machos p 



I 



^-TTU 2 T Rl{s)ds. (4) 



m 



Unlike the rate or average duration, it is independent of the Macho mass distri- 
bution. For this reason it also contains no direct information about the mass of 
the lensed objects. Note also that since detection efficiencies depend upon the 
duration of events, it is important to have models which predict durations. Now 
to find the rate at which Machos enter the microlensing tube requires knowledge 



of the distribution of velocities all along the tube. So the rate and the distri- 
bution of event durations are hard to calculate because they require the entire 
phase space distribution function (DF) i^(v,x). The differential rate is given by: 

1 3 

dT = — F(v, x) cos 0UTR e v±_d vdxda, (5) 

m 

where the angles and notation are defined in Griest (1991). 

What makes the calculation particularly difficult is that the DF cannot be 
prescribed arbitrarily as a Maxwellian, for instance. This is because the Machos 
are collisionless, so the DF is constrained to obey the collisionless Boltzmann 
equation. By Jeans' theorem, this implies that the DF depends only on the iso- 
lating integrals of motion (see Binney & Tremaine 1987, p. 220). Self-consistent 
solutions for distributions of velocities that build flattened halo models are scarce. 
The largest known set of axisymmetric models with simple DFs are the "power- 
law galaxies" (E93, E94). These form the basis for the exploration of microlensing 
in this paper, allowing us to go beyond simple spherical models. Note, however, 
that all these models are axisymmetric and oblate, while N-body simulations 
suggest that halos may well be triaxial. Exploration of triaxial models will be 
done in a future paper. 

The parameters of the power-law models are: 

(1) The core radius R c , which measures the scale at which the density law begins 
to soften. 

(2) The flattening parameter q, which is the axis ratio of the concentric equipo- 
tential spheroids, with q = 1 representing a spherical (E0) halo and q ~ 0.7 
representing an ellipticity of about E6. The "isophotal" ellipticity of the dark 
halo is a function of q, as well as other parameters of the model [see E94, eq. 
(2.9)]. 

(3) The parameter (3, which determines whether the rotation curve asymptotically 
rises, falls or is flat. At large distances R in the equatorial plane, the rotation 
velocity v c [ rc ~ R~ ' . So j3 = corresponds to a flat rotation curve, while j3 < 
is a rising rotation curve and fi > is falling. 

(4)The solar radius i?o, which is the distance of the Sun from the galactic center. 

(5) Finally, the normalization velocity vo, which determines the overall depth 
of the potential well and hence the typical velocities of Machos in the halo. In 
the limit fi = 0, q = 1 and large R (spherical halo with a flat rotation curve) 

"0 = ^circ- 

Using z as the height above the equatorial plane, the potential of the power- 
law models is 



v^c/(3 



*= I 



{Rl+R 2 +z 2 q- 2 )P/ 2 ' 



,2 



if /MO, 

(6) 



.^-lo g (R 2 c+ R 2 + z 2 q- 2 ), if = 0, 



(7) 



and the mass density is 

_ vlRP R 2 C (1 + 2g 2 ) + R 2 (l - (3g 2 ) + z 2 (2 - (1 + (3)g- 2 ) 
9 4jrGq 2 (Rl+R 2 +z 2 q- 2 )(P+ 4 )/ 2 

The DF corresponding to this potential-density pair is 

' A^|£|4//3-3/2 + £|£|4//3-l/2 + c ^2//3-l/2^ if /? ^ 0, 

F(E,L Z )= < 

, AL 2 z exp(4:E/v%) + Bexp(4:E/v%) + Cexp(2E/v%), if (3 = 0. 

where the constants A, B, and C are given in E93 and E94. As required by Jeans' 
theorem, the DFs depend only on the isolating integrals of motion, namely the 
relative energy per unit mass E = \P — ^v , and the angular momentum per unit 
mass about the symmetry axis L z . The circular velocity in the equatorial plane 
is 

„,2 d/3 e>2 



»L = ,„^tL, w - P» 



(B- + iJ 2 )W+ 2 )' 2 ' 

Note that the limit q = 1, fi = and R c = recovers the standard singular 
isothermal sphere used by Paczyhski. Allowing a core radius gives 

Jltf+Mg (10) 

This differs from the cored isothermal sphere considered by Griest (1991) in sev- 
eral ways. First, the rotation curve approaches its asymptotic value more quickly. 
Second, the DF given by the q = 1, fi = limit of equation (8) is self-consistent, 
whereas Griest (1991) assumed an approximate Maxwellian distribution of ve- 
locities. 

So far we have only modeled the dark halo. However, in the standard model, 
a substantial fraction (~ 40%) of the centripetal force at the solar radius derives 
from the disk stars. This is represented by a thin exponential disk with a scale 
length of Rfj = 3.5 kpc, normalized to a surface density of So = 5OM0pc at the 
solar radius (Gilmore et al. 1989; Gould 1990). It is possible that the disk of our 
galaxy is substantially larger than the canonical value. (Oort 1960; Bahcall 1984; 
Kuijken and Gilmore 1989; Gould 1990). Recent microlensing results (Alcock, 
et al. 1994; Udalski, et al. 1994) as well as studies of the optical rotation curves 
of external galaxies (Buchhorn 1992; Kent 1992) may suggest this. We consider 
such a "maximal disk" by taking So = IOOMqpc . The rotation velocity added 
in quadrature is thus (Freeman 1970; Binney & Tremaine 1987, p. 77) 



disk 



AnG^hy 2 [I {y)K {y) - h{y)Ki{y)] , (11) 



where y = R/(2Rd), and the /„ and K„ are modified Bessel functions. Note 
that in adding a contribution from the disk to the local circular velocity, we have 



sacrificed self-consistency. Really, we should find the DF of the power-law halo 
in the combined potential field of both disk and halo - instead, we use the DF 
(8). As has been argued elsewhere (Evans & Jijina, 1994), this is an reasonable 
approximation for the LMC, SMC and M31, where microlensing typically occurs 
at heights above the equatorial plane of many kpc. 

In this paper, our aim is to estimate the contribution of the Galactic halo to 
microlensing. Of course, this is not the only possible source of deflectors. Towards 
the LMC, there is the possibility of microlensing by the LMC dark halo or disk 
(Gould 1993b, Sahu 1994). The optical depth is ~ 2.5 x 10 -7 for microlensing by 

LMC halo lenses and ~ 0.09 X 10 for LMC disk lenses. The Galactic halo makes 
a contribution that is roughly three times greater and so is the dominant source 
of lenses. However, this is not the case for lines of sight towards M31. The optical 
depth is dominated by Machos in the halo and disk of M31 (Crotts 1992, Gould 
1993a). Crotts (1992) estimates that the halo of our own Galaxy contributes just 
20% to the total optical depth. Microlensing towards the Galactic bulge poses 
perhaps the hardest problems of separating the contributions of different deflector 
populations. Bulge stars can undergo microlensing not only by halo Machos, but 
also by other bulge and disk stars (Griest et al. 1991, Paczyhski 1991, Kiraga 
& Paczyhski 1994). At Baade's Window, the optical depth is ~ 6.3 X 10~ 7 for 
microlensing by bulge lenses, ~ 5.0 X 10 for disk lenses. The dark halo only 
makes an important contribution if the core radius is small. 

We are now in a position to calculate the microlensing observables - the 
optical depth, rate and average duration of events. They can be found using 
equations 2, 3, 4 and 7. The results are single quadratures and readily evaluated 
on the computer. They are displayed in Appendix A. In Appendix B, we give 

the differential microlensing rate dT/dt, where t is defined just after equation (1). 

The probability of obtaining an event of duration t is just (dT/dt)/T. 

3. Range of models 

In order to explore the scatter in microlensing observables, we build a set 
of halo models which span the observationally allowed range. The power-law 
galaxy models allow us to vary the flattening, core radius and rotation law, and 
we consider both canonical and maximal disks. For each parameter in the model, 
we therefore find the range permitted by the observations. Then, several values 
of each parameter are chosen to represent the range. We also ensure that each set 
of parameters gives a model consistent with the measured Milky Way rotation 
curve. So, we study the statistical properties of an ensemble of models, each one 
of which is a plausible representation of the dark halo of the Milky Way. 

For the dark halo flattening, little is known. So the entire range of flatten- 
ing allowed by the power-law models is examined. This varies between E0 or 
spherical (q = 1) and roughly E6 or E7 (depending on (3). The core radius of the 
dark halo is also uncertain - Bahcall, Schmidt & Soneria (1983) estimate R c as 
2 kpc from star count data, while Caldwell & Ostriker (1981) suggest 10 kpc. If 
the disk is maximal, values as large as 20 kpc are possible. We consider values 
of 2 kpc, 5 kpc, 10 kpc, and 20 kpc. The parameter fi determines the slope of 
asymptotic circular velocity. Between Rq and 2i?o, the circular velocity is prob- 
ably within 10 — 15% of the I.A.U value of 220 km/s, but whether the measured 
HI rotational velocities rise or fall with R depends upon estimates of the solar 
position Rq and the local circular speed v circ (Ro) (see Fich et al. 1989; Jones et 
al. 1993). Beyond 20 kpc, little is known directly, though arguments based on 
the kinematics of distant satellite galaxies support the idea of a relatively flat 



rotation curve out to an unknown cut-off (Fich & Tremaine 1991). However, 
current theories of galaxy formation tend to favor the alternative view that dark 
halos extend indefinitely, fading into structure on larger scales. So, we do not 
consider a halo cut-off in this paper - it would add yet another poorly known 
parameter to our model. We investigate power-law halos with fi = -0.2, 0, and 
0.2. These correspond to rotation curves which rise by ~15%, are flat, or fall 
by ~15% between the solar radius and twice the solar radius, depending a little 
upon R c . 

The value of the solar radius Ro has been reviewed by Reid (1989). He 
shows that most recent determinations lie between 7 kpc and 9 kpc, with 7.7 
kpc being his preferred value. This differs considerably from the IAU value of 
8.5 kpc (Kerr & Lynden-Bell 1986). We examine the values Ro = 7, 8, and 9 
kpc. Finally, perhaps the single most important parameter is the normalization 
velocity vq. Given our fixed disk contribution to the total rotation law, the 
parameter vq is now specified once we settle upon a choice for v c i rc (Ro). Merrifield 
(1992) estimates v circ (Ro) = 200 ± 10 km/s, Fich, Blitz, & Stark (1989) give 
Ucirc(-Ro) = 220 ± 30 km/s, while Rohlfs et al. (1986) give values between 170 
km/s and 200 km/s between Rq = 6 kpc and Rq = 16 kpc. For our ensemble of 
models, we impose the constraint that the total circular velocity lies between 180 
km/s and 250 km/s at Rq and 2Rq. Note that the IAU value is 220 km/s (Kerr 
& Lynden-Bell 1986). We also investigated a more restricted ensemble of models 
with 190 < v c i rc (Ro) < 230 km/s. We find all our results also hold for this more 
restricted ensemble. 

4. Uncertainties in the Rates 

First, let us consider the difference caused by using the optical depth instead 
of the microlensing event rate. The optical depth to microlensing is the mean 
number of Machos in the microlensing tube; that is the number of microlensing 
events taking place at a given moment. It is easy to calculate since it is inde- 
pendent of lens mass and velocity, and only requires knowledge of the density 
distribution /3(x). For this reason, it is the most widely estimated quantity. But 
how well does it trace the microlensing rate? 

We are able to answer this question since both the rate T (equation Al) and 
the optical depth r (equation A6) are known for the power-law models. One way 
to test this is to plot (7), which is the ratio of optical depth r and T, ft) = —t/T, 
for many different models. The average duration (tj is a constant if r and T are 
well-correlated. In Fig. 1, we show histograms of (7) for microlensing towards 
the LMC, SMC and M31 for our ensemble of models. Figs, la-c demonstrate 
that (tj tends to vary by more than a factor of two between models. Figs, ld-f 
show an even larger spread for maximal disk models. Figs. 2a-f show this another 
way by plotting the rate vs the optical depth for the set of models. These plots 

show that ft) is indeed much less model dependent than either r or T. While the 
rate and r vary by more than a factor of ten in these plots, their ratio varies only 
~ 2 for a canonical disk. In fact, we note that the line T oc t ' is a fairly good 
fit to all the models we have considered. Thus the large scatter in ft) seen in 
the maximal disk histograms is mostly just due to the large scatter in rate. (The 



~k To the extent that the relation T = ar) m f 1 ' 2 r 3 ' 2 holds, where a is a constant from theory, 

we have that a = / _1 ' 2 77~ 1 (4/7r) 3 ' 2 (t) r -1 ' 2 is independent of the model parameters. 

Thus, if the macho fraction / were known, one could extract the mass integral r\ m from 



rate varies more than the optical depth.) In all the plots we use m = Mq, but for 
an arbitrary mass distribution just scale T by eta m and ft) by ry~ . Keep in mind 
that a given experiment can produce only one point in the T, r plane and that 
the primary use of a measurement will be to find /, the Macho fraction. We see 
that for approximate work, the optical depth does a reasonable job of predicting 
the rate. But for more detailed work, especially when efficiencies are involved, 
the difference between rate and optical depth should be kept in mind. We also 
note that the predicted distribution of event durations is found as a differential 
rate (Appendix B). 

Next let us turn to scatter in the predicted microlensing rate caused by un- 
certainties in the halo parameters. Fig. 2 shows that for all lines of sights, there 
is a scatter in the rate of more than a factor of ten for a canonical disk. For 
the LMC, the models with the smallest rate have spherical halos with small core 
radii, falling rotation curves, and small values of vo, while the models with the 
largest rates have either spherical or flattened halos, but large core radii, rising 
rotation curves, and large values of vq. This is as expected, since any model 
which puts more mass at a large distance in the direction of the LMC will have 
a larger microlensing rate, and a larger optical depth. This is shown in Fig. 3 
in which we plot the optical depth against the rotation velocity at r = 50 kpc. 
The correlation between r and v c {j- = 50 kpc), while not perfect, is quite good. 
Note that the mean value of the rotation velocity at Ro is nearly independent of 
the microlensing rate. Figs. 2d-f and Fig. 3c-d show the case of a maximal disk. 
Here we see that the rate and optical depth can be considerably smaller than 
for a canonical disk. Also there is a variation between models of several orders 
of magnitude. This is as expected since in these models the disk is the main 
contributor to the rotation curve at the solar distance. Thus a smaller enclosed 
halo mass is required to match observations, and the halo parameters are poorly 
constrained. 

The halo may only consist of a fraction / of baryonic matter in the form 
of Machos. Thus, a factor of more than ten uncertainty in the predicted rate 
caused by the poorly determined halo parameters makes it difficult to determine 
the allowed amount of non-baryonic dark matter. It is clearly essential to reduce 
the uncertainty. 

5. Reducing Model Uncertainties 

The primary way of reducing the model uncertainties in the microlensing 
observables is to determine the halo parameters. Even within the restricted 
framework of the power-law galaxy models, if (3, vq, q, Ro and R c are known, 
there is still uncertainty in the rate. This is because the DFs equation (8) are 
the simplest consistent with the potential and the density, but are certainly not 
unique. There are still further DFs that depend on non-classical third integrals 
of motion and generate anisotropic velocity distributions. Note, too, that even 
though our models give a plausible representation of the Milky Way, there cer- 
tainly exist other alternatives (see e.g., Frieman & Scoccimarro 1994; Gates & 
Turner 1993; Giudice, Mollerach & Roulet 1994) with different lensing properties. 
And of course, the size of the disk plays a crucial role. 



observables r\ m k, / 1 / 2 a~ 1 (4/7r) 3/2 # e //£ ,1/2 (X^'/£i)~ 3/2 > where E is the total exposure, 
N e ff = Y^ e 7 ! e i i s the efficiency at which events of duration t{ are recovered, and the sums 
go from 1 to the number of observed microlensing events. For LMC microlensing in our set 
of models we find a « 3850 ± 260yr -1 . The physical basis for this relationship may simply 
be that the optical depth is proportional to the mass along the line of sight ex i> 2 , and the 
rate is proportional to the optical depth times v c . 



One obvious way to determine halo parameters is to use conventional astro- 
nomical techniques - observations of stars, gas and satellites - to fix the solar 
radius and circular speed more accurately. For example, fixing the solar radius 
at 8 kpc, and demanding v c { rc = 220 km/s ±5% between 8 and 16 kpc reduces 
the spread in microlensing rates toward the LMC from more than a factor of ten 
to a little more than a factor of two (for the canonical disk). Uncertainties in r 

and ft) are reduced similarly. A better determination of the halo core radius by 
stellar observations would also be important. 

However, it is also possible to use the microlensing experiments themselves to 
determine the halo parameters and reduce the model uncertainty. The basic idea 
is to exploit the fact that there are at least four viable lines-of-sight out of the 
Milky Way in which to measure the microlensing rate and average event duration. 
Each line-of-sight (LMC, SMC, M31 and the bulge) offers a different "pencil 
beam" through the dark halo, and so by comparing the rates, optical depths, 
and average durations among the different lines-of-sight information concerning 
the halo shape can be gained. Several of the parameters, such as flattening q 
and asymptotic slope of the rotation law (3, may best be determined this way. 
So microlensing gives us a new probe of the density and velocity structure of the 
dark halo. This is in addition to information on the size of the disk gained via 
microlensing. 

For instance, a scatter plot of the ratio of LMC and SMC rates vs the LMC 
and SMC average durations is shown in Fig. 4a. The models clearly fall into 
two distinct groups. Those models marked with a circle all have round halos 
(E0), while those with a square are flattened to roughly E6. Thus the ratio of 
LMC rate to SMC rate is a excellent indicator of halo flattening. This effect was 
first discovered - using optical depth rather than microlensing rate - by Sackett 
& Gould (1993). Frieman & Scoccimarro (1994) have recently cautioned that 
the robustness of this diagnostic may be lost if the halo is tilted with respect to 
the disc - although such a configuration cannot be a long-lasting equilibrium. 
So, the halo flattening can probably be determined if enough events are found 
to allow accurate measurement of the SMC microlensing rate. Figs. 4b and 4c 
show the rate ratio for M31/LMC, and M31/SMC. While separation of flattened 
models is still evident, one sees from the figures that it is the LMC and SMC 
position relative to the halo axis of symmetry that make the measurement of 
the flattening so easy. Note again, that in an experiment one measures only one 
LMC rate (and optical depth) and one SMC rate, and so gets only one point in 
any of these scatter plots. It is also interesting to observe from Fig. 4a that the 
model uncertainties in the LMC/SMC rate ratio are much greater for flattened 
halos than for spherical halos. The case of a maximal disk is not shown, since it 
looks almost identical to Fig. 4. 

Can we use microlensing to determine whether the halo has a rising or falling 
rotation curve? The LMC and SMC are at nearly the same distances (50 and 60 
kpc), so it is natural to expect the ratio of M31 to LMC microlensing to be the 
most useful discriminant. Note that rate ratios are convenient to use, because the 
magnitude of any rate always contains the unknown parameter /. In Fig. 5 we 
plot the M31/LMC rate ratio vs the LMC rate for the set of models above, with 
triangles for (3 = 0.2 (falling rotation curve), circles for (3 = (asymptotically 
flat rotation curves), and stars for (3 = —0.2 (rising rotation curve). In Fig. 5a, 
all models are plotted, while in Fig. 5b and Fig. 5c only models with spherical 
(q = 1) and flattened (q = 0.71 or q = 0.78) halos respectively are shown. 
In Fig. 5a some separation of models with different values of (3 is evident but 
there is substantial ambiguity, which would make a direct estimate of (3 using 
this method difficult. However, suppose that we have already determined the 
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halo flattening by use of the ratio r LMC /r 'sMC- Then, as shown in Figs. 5b 
and 5c for a canonical disk, a fairly clear separation of rising, falling, and flat 
rotation curve parameter can be accomplished. Thus, the ambiguity seen in 
Fig. 5a is largely removed when models with different flattenings are plotted 
separately. The exception is some overlap between models with R c = 2 kpc and 
R c = 20 kpc and different values of fi. This ambiguity is probably removable 
as discussed below. The case of a maximal disk is not displayed, as it is very 
similar. So, the asymptotic form of the rotation law, or equivalently (3, can 
probably be determined from the M31/LMC rate ratio once q is known. Keep in 
mind, however, the caveats mentioned in §3 concerning our M31 rate calculation, 
which may result in corrections which modify this effect. If halo microlensing 
can be distinguished from M31 microlensing, a measurement of fi should then be 
possible. 

Next, can we determine the halo core radius R c ? The parameter R c affects 
mainly the inner portion of the halo and overall normalization of the halo mass. 
This overall normalization is mixed in with vq and /, and so the best hope 
in determining R c is probably a comparison of the bulge with a more distance 
source such as the LMC. Here we have the problems mentioned in §3 concerning 
bulge microlensing; our modelling of the distribution of velocities is not adequate 
along the disk. But, the optical depth is independent of the velocities and will 
give some indication of the rate. Even so, our calculations do not give the total 
optical depth towards the bulge, merely the contribution of the optical depth 
from the halo. 

In Fig. 6, we plot the LMC/bulge optical depth ratio vs the bulge optical 
depth, where triangles indicate R c = 2 kpc, boxes indicate R c = 5 kpc, and 
stars indicate R c = 10 kpc. A reasonably clean separation is obtained when 
this ratio is plotted for all the models (not shown). In Fig. 6, this separation 
is made clear-cut, if one supposes fi and q have already been measured by the 
methods above. The R c = 20 kpc models have an LMC/bulge ratio of greater 
than 10, and are very easily distinguished even with no prior knowledge of fi 
and q. (They fall off the top of the plots in Fig. 6). Even if fi and q are not 
known, the separation is quite good if the value of the solar radius Ro is held 
fixed. So, a better determination of Ro by non-microlensing means can allow a 
clearer separation of the effect of the halo core radius. The case of a maximal 
disk is not shown since it gives very similar results. 

6. Distribution of Event Durations 

Since the duration of a microlensing event is proportional to the Einstein 
radius (oc m ' ), the duration of an event gives information about the mass of lens 
which caused it. In trying to understand the nature of the objects responsible for 
the observed microlensing, this is important information. But the duration also 
depends upon the unknown lens velocity and distance. Thus, a given mass Macho 
can cause a wide distribution of event durations. This distribution must be used 
statistically to infer probable masses from observed durations. Using the DF's 
(equation (8)), the distribution of event durations can be found. The formula and 

definitions are given in Appendix B. In Fig. 7, we show several t distributions. 
One sees that different halo parameters give quite different distributions. It is the 
average of these distributions ft) that is shown in the histograms in Fig. 1. Fig. 7 
shows that, as expected, uncertainty in the halo model will lead to additional 
uncertainty in determining the masses of the lensing objects. The curves labeled 
(a), (b), and (c) are canonical disk cases with various choices of halo parameters, 
while curve (d) shows a maximal disk example. We also note that the scaling 
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introduced in Griest (1991) works fairly well for models we considered. That 

is, by scaling the t axis by ft) , and the dT/dt axis by (t\ all the curves are 
found to lie roughly on top of each other. This means that for power law galaxy 
models along a given line- of- sight, the shape of the distribution is much more 
model independent than peak value. 

In a future paper we plan to explore further the information that can be 
extracted from event duration distributions, and include other possibilities such 
triaxiality, streaming motion, etc. 

7. Conclusions 

This paper has shown how to exploit the power-law galaxy models (E93, E94) 
as simple, flexible and realistic representations of the dark halo. These models 
have the advantage of simple and analytic phase space distribution functions and 
therefore permit accurate calculation of the optical depth, microlensing rate and 
average event duration. We provide formulae for these quantities as a function of 
the halo parameters and source distance and direction (Appendix A). The distri- 
bution of event timescales is presented in Appendix B. We apply our formulae to 
study microlensing towards the Large and Small Magellanic Clouds (LMC and 
SMC), the galactic bulge, and the M31 disk galaxy. We find that: 

(1) For a canonical disk, the optical depth is a reasonable indicator of the mi- 
crolensing rate to within a factor of two. This is important, because the optical 
depth is much easier to calculate than the rate and probably will continue to be 
widely used by investigators. For more accurate work, as well as for derivations 
of the distribution of durations, galaxy modeling with distribution functions is 
crucial. For a maximal disk the agreement between optical depth and rate is less 

robust, though the relation T oc m ' t ' seems to hold. 

(2) The evaluation of the fraction / of the halo consisting of Machos is hampered 
by the uncertainties in the galactic constants, such as the shape of the rotation 
law and the flattening of the dark halo. For a realistic set of halo models, we 
found rates toward the LMC and SMC can vary by more than a factor of ten 
from model to model for a canonical disk, and by several orders of magnitude for 
a maximal disk. Left unaddressed, this model uncertainty will thwart accurate 
determination of /. 

(3) An attractive way of reducing the uncertainty - which simultaneously opens 
up a new method in galactic astronomy - is to use microlensing to explore the 
shape and structure of the dark halo. This has also been realised by Sackett & 
Gould (1993), who showed that the ratios of optical depth towards the LMC and 
SMC is a robust indicator of the flattening of the dark halo. We confirm this 
result by showing that the ratios of the event rates also distinguish flatness. In 
particular, the ratio of microlensing rates towards the LMC and SMC is ~ 0.7—0.8 
for E0 halos and ~ 1.0 — 1.2 for E7 halos. This is true for both canonical and 
maximal disk models. Once the flattening has been established, the asymptotic 
slope of the rotation curve fi might be determined using the M31/LMC rate ratio. 
The LMC/bulge ratio contains important information on the halo core radius. 
We caution that this may not be easy to extract, as the dark halo is probably 
not the dominant source of lenses towards the bulge. 

In summary, the discovery of a dark halo consisting of a significant fraction 
of Machos is only the starting point for an exploration of the halo characteristics 
which microlensing can help determine. 
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Appendix A 

In this appendix, we give the formulae for the microlensing rate and optical 
depth for the general flattened halo model described in the text (equations 4-7). 
The total rate T of microlensing in a power-law halo with model parameters 
f3 } v 0} R C} R , and q is 

Cqu t t,giff /2 (/? + 2)(l-g 2 )rM 

y/2irM/M @ 2« z 2 ^£ 1 / 2+3/3 / 2 r(^) r 



Cqu t v 3 R 2 c + ^ /2 ((3 + 2)T(n p ) 
^ttM/M© cq 2 ^^! 2 ^! 2 T (dp) 2 

Cou T v 3 Rf /2 (2- 1 -^) T (np-2/\f3\) 

\/27tM/M q c^Li/2+3/3/2 T(d f3 -2/\(3\) 3 ' 

Here, Co = 1/s/GMq, T(x) is the gamma function, and the integrals I{ are 
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. dsy/s{l - s)(A's 2 + B's + C) 



I 



{D's 2 + Es' + F'f+WI 2 
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f dsy/s{l- S ) 

J (D's 2 + Es' + F') 2 +^/ 4 { ' 



o 
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J (D'i 



dsy/s(l — s) 



with 



(D's 2 + Es' + F') l +^/ 4 ' 
o 



A =3 cos 6, B = — 6Ro cos bcos£/L, 

C' =2Rl/L 2 + Rl cos 2 £/L 2 + R 2 sin 2 £ sin 2 b/L 1 

D = cos b + q~ sin 6, E = — 2Rq cos bcos£/L 

F'=(R 2 C +R 2 )/L 2 . 



■2 



(A3) 



The quantities 6, / are the galactic coordinates of the source star, L is the source 
distance, G is Newton's constant, and c is the speed of light. The constants np 
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and dp have a different form according to whether fi is positive or negative 



n P = \ 



-4 3 



/? 



U 



if (3 < 0, 



(A4) 



+ 2, if (3 > 0, 



da = < 



— -1, if/?<0, 

P 



(A5) 



In the limit /? — )■ (the case of an asymptotically flat rotation curve), the ex- 
pression for the rate follows from the above by systematic use of the formula 

T(x + l/2)/r(x) — )■ ■sjx asi-y oo. The optical depth r is 



,2z?/3 2 



v^BicU l T f s(l-s)(A"s 2 + B"s + C")ds 



I 



c 2 q 2 L? J (D's 2 +E's + F l )^+ 4 )/ 2 



(A6) 



where 



A" = (1 - (3q 2 ) cos 2 6 + (2 - (1 + (3)q~ 2 ) sin 2 6, 

5" = -2(1 - f3q 2 )R cos 6 cos ^/L, 

C" = (R 2 c (l+2q 2 ) + R 2 (l-(3q 2 ))/L 2 . 



(A7) 



The quadratures are straightforward to evaluate on the computer. 



Appendix B 



The distribution of event durations is important for finding the mass of 
the lensing objects. It is given by the normalized differential microlensing rate 

(dT/dt)/T, where t = 2R e /v_\_, and v±_ is the speed of the Macho perpendicu- 
lar to the line- of- sight. The time the Macho spends inside the Einstein radius, 
t e = {uPp — u mi ^] I ' t, where u mm is defined in equation (1), and uy = 1. The 
average duration is related to the average t by (t e ) = j (t). In many cases 

it is advantageous to use distributions in t, since they are independent of the 
amplifications. 
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For the model described in the text, we find: 



dT _ u T f L 6 
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and the integrals J 8 - are 
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If fi < 0, the integrals are evaluated over the interval [0, 1]. If fi > 0, then we 
must restrict the domain of integration by 



- 2 8(iL s(l - s) (m/M®) (D 1 s 2 + E's + F'fl 2 



(v Q cCo) 2 (Rc/Lf 



(B4) 



The constants A', B' , C", -D' and E' are given in Appendix A. The additional 
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constants are 
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Figure Captions 

Figure 1: Histograms of the average duration (tj = —t/T for the ensemble of 
halo models discussed in the text. Part (a) is for the LMC, (b) is for the SMC, 
and (c) is for M31. If optical depth tracked microlensing rate perfectly each 
histogram would be a delta function. Parts (d)-(f) are the same for a maximal 
disk model. Note all plots are for m = 1M©; scale by ry~ for other masses 
(Equation (3)). 

Figure 2: Scatter plots of microlensing rate vs optical depth for the ensemble 
of models discussed in the text. Part (a) is for the LMC, (b) is for the SMC, and 
(c) is for M31. Each point represent a consistent model of the dark halo. Parts 
(d)-(f) are the same for a maximal disk model. All event rates scale Y oc r/ m . 

Figure 3: Scatter plots of optical depth vs v^(50 kpc), the square of the total 
rotation velocity at 50 kpc in the galactic plane. The mass of the Galaxy interior 
to this distance is proportional to this squared velocity. Parts (a) and (b) are for 
a canonical disk, while parts (c) and (d) are for a maximal disk. 

Figure 4: Finding the flattening parameter q. Scatter plots of the ratio of 
rates vs. the ratio of event durations. The circles represent halo models which 
are spherical (q = 1), while the squares represent flattened halos (q = .71 for 
(3 = 0, -0.1; q = 0.78 for (3 = 0.1). part (a) is for LMC/SMC and shows clear 
separation of spherical and flattened halos. Part (b) is M31/LMC, and part (c) 
is M31/SMC. 

Figure 5: Finding the asymptotic slope (3. Scatter plots of the M31 rate divided 
by the LMC rate vs the LMC rate. The stars represent halo models with (3 = —0.2 
(rising rotation curve), the circle models with (3 = (flat), and the triangles 
models with (3 = 0.2 (falling). Part (a) shows all models, while part (b) shows 
only spherical models and part (c) only the flattened models. Separation of the 
models becomes easier if the flattening is known. The line of ambiguity in some 
panels is due to R c = 20 kpc models, which can be distinguished as shown in 
Figure 6. All event rates scale T oc r/ m . 

Figure 6: Finding the core radius R c . Each panel shows models of definite values 
of (3 (rotation curve slope) and q (flattening). The separation between models is 
quite good. The triangles represent R c = 2 kpc, the squares represent R c = 5 
kpc, and the stars represent R c = 10 kpc. Panels marked E0 are for spherical 
halos, while those marked E6 are for flattened halos. Models with R c = 20 kpc 
were also considered but they are easily distinguished since they typically have 
r(LMC)/r(bul) > 10 and therefore fall off the top of the figures. 

Figure 7: Examples of LMC t distributions for various model parameters. The 
integral under each distribution is unity. The curve marked (a) is for a "standard" 
spherical halo ((3 = 0, q = 1, R c = 5 kpc, vq = 200 km/sec, Rq = 8.5 kpc, and 
r = 1.64 X 10 events/yr). Curve (b) has a shorter average duration ((3 = —0.2, 
q = 1, R c = 5 kpc, v Q = 200 km/sec, R Q = 8.5 kpc, and T = 3.9 x 10 -6 
events/yr). Curve (c) has a longer average duration ((3 = 0.2, q = 0.78, R c = 10 
kpc, vq = 210 km/sec, Rq = 8.5 kpc, and T = 1.24 X 10~ 6 events/yr). Finally 
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curve (d) has a maximal disk, which greatly reduces the amount of halo material 
(/? = 0, q = 1, R c = 20 kpc, v Q = 90 km/sec, R Q = 7 kpc, and T = 9.37 x 10 -8 
events/yr). The average of each distribution is (7). All event rates scale T oc r] m . 
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